A note on a relation between ac Josephson effect and double-well BEC oscillations 



o 
o 

(N 
O 

Q 



C/2 ■ 

ctf ■ 

^— > ■ 
C ■ 
c3 ■ 
=5 

+-» ■ 



i 

c 

o 
o 



> 

00 

m 
m 

o 



X 



Leo Radzihovsky and Victor Gurarie 
Department of Physics, University of Colorado, Boulder, CO 80309 
(Dated: December 17, 2009) 

In this brief note we comment on the relation between the ac Josephson effect and the coherent 
oscillations of a Bose-Einstein condensate confined to a double- well potential. The goal is to elucidate 
the extent to which the latter is a realization of the former. We detail the correspondence that 
emerges in the high occupation limit of the double-well potential, and particularly note the relation 
between the two oscillation frequencies. 



A realization of trapped degenerate atomic gases has 
opened opportunities to study many interesting quan- 
tum many-body phenomena in previously unexplored 
regimes lj. Recent experiments @, Q on Bose-Einstein 
condensates (BEC), trapped and oscillating in an un- 
balanced double-well potential have sought to realize an 
atomic Bose gas analog 0, H, H, 0, [1] of the alternating- 
and direct-current (ac and dc) Josephson effects [9l. lioj]. 

While there are some obvious analogies, considerable 
fundamental differences between the two systems exist. 
For instance, a conventional Josephson junction ( JJ) be- 
tween two superconductors is an open quantum many- 
body system driven by a fixed electro-chemical poten- 
tial (voltage) difference between the left and right con- 
tacts, or with a current imposed by an electrical circuit. 
Thus on general grounds (gauge invariance and Hciscn- 
berg equation of motion) the evolution of the relative 
phase <j) — 4>l — 4>r in an ac Josephson effect is given, 
exactly, by 



hcf) = 2eV, 



(1) 



where eV = [II — ^r is the imposed electro-chemical 
potential difference across the left and right contacts 
controlled by the voltage V. When combined with the 
(lowest harmonic) expression for the Josephson current, 
I = Jo sin <^), the above exact linear growth of <p(t) with 
time gives the standard ac Josephson effect, 



lit) = 7 S m(ujjt), ujj = 2eV/H. 



(2) 



In addition to its basic physical importance it pro- 
vides an extremely accurate voltage-frequency relation 
(thereby defining a Volt to one part in 10 s ) with a num- 
ber of other important applications[10j. In a conven- 
tional superconducting JJ the critical current is given 
by the standard Ambegaokar-Baratoff formula 11], Io — 
7rAG„/(2e), where A is the superconducting (pairing) 
gap and G n is the normal state Josephson junction con- 
ductance. 

In contrast, a double- well BEC is a closed system with 
only the total number of atoms fixed and is not necessar- 
ily in the thermodynamic limit. Deep in the condensed 
state and for deep wells the coherent dynamics can be 
studied via two coupled (Gross-Petaevskii) equations of 
motion for two coherent-state amplitudes 5 , L,fl(i), with 
|^l,_r| 2 = NL,R.{t) giving the number of atoms in the left 
and right wells, with only the total number N — Nl+Nr 



of atoms conserved. The dynamics is given by Euler- 
Lagrange equations for the coherent-state action given 
by 



S 



dt [** L ihd t *L + y R ihd t y R - H(? L , , (3) 



with the Hamiltonian 



e L \* L \ 2 +e R \V R \ 2 (4) 



9o_ 

2JV' 



go 

2N' 



and parameters J,£L.R-,9i straightforwardly derivable 
from a continuum model of interacting bosons trapped in 
a double- well potential Since the Hamiltonian should 
scale linearly with the system size, for convenience we 
defined the interaction couplings with explicit factors of 
1 /N so that the parameters gi do not scale with the sys- 
tem size. 

For gi = the dynamics reduces to that of a two-level 
system (e.g., spin in a magnetic field), with oscillations 
(that we will loosely call Rabi oscillations) arising from 
non-eigenstate initial conditions. This two level system 
has a matrix representation 



H 



sl -J 
-J £R 



with well-known eigenvalues 



E± = e ± \/j 2 + h 2 , 



where we defined 



£r), h = -{e L 



sr) 



(5) 



(6) 



(7) 



(with the energy difference h not to be confused with the 
Planck constant, for which we will use the symbol 2ttH 
throughout). Thus any quantity quadratic in ^ l and ^ R 
will oscillate with the Rabi frequency 



ujr = (E + - E_)/h = 2y/j 2 + h 2 /h. 



(8) 



Although this double-well BEC system was studied in 
a seminal work by Smerzi, et al.Q and a number of works 
that followed, from the recent discussion [l2j of the lat- 
est experiment Q it appears that a number of questions 



2 



remain unanswered, namely: (i) How is the Rabi fre- 
quency, lor, which depends on system-specific quantities 
such as J, related to the universal Josephson frequency, 
uij in which depends only on the applied voltage? (ii) 
Moreover, what role do interactions (clearly neglected in 
ljr) and the thermodynamic limit play in establishing 
the relation between the two systems and corresponding 
frequencies? 

A connection between these (otherwise quite distinct) 
systems only exists in the specific limit of macroscopic 
wells, Ll,r 00 (reached for large occupation N L>R ^> 
1) and a finite barrier thickness d -C L of the double- well 
BEC system. This is necessary in order to approximately 
model the thermodynamically large lead reservoirs of a 
superconducting Josephson junction. Our key observa- 
tion is that in this limit the Josephson coupling, J, as 
defined by ([4} vanishes as l/L, i.e., vanishes in the large 
atom occupation number limit, and thus ljr goes over to 
u)j, exactly in the thermodynamic limit, with the iden- 
tification of h with eV. Moreover, in this limit the am- 
plitude of oscillations in the number imbalance is always 
small regardless of the interaction strength, and thus the 
nonlinearity plays negligible role in their dynamics, aside 
from just redefining the relationship between h and V to 
be 



eV <-> h+ {g Q - gi) 



N 



(0) 



N 



(0) 



2N 



(9) 



with N^ R the equilibrium number of atom in the two 
wells. 

The vanishing of J for macroscopic leads and fixed 
barrier width d can be understood on general grounds 
simply by noting that J is an interfacial energy per par- 
ticle, associated with the coupling of the left and right 
leads. It is therefore proportional to the surface area 
of the barrier A ~ L 2 divided by system's volume, i.e., 
J ~ L 2 £/L i ~ £/L, (with £ the inverse penetration length 
scale set by the barrier width and height) and thus indeed 
vanishes in the above thermodynamic limit. 

To derive these results more explicitly, we estimate 
the typical size of J in a macroscopic system by solv- 
ing the Schrodinger equation of a particle of mass m a in 
a double-well symmetric potential, U(z) — U(—z), with 
the difference between the first excited and ground states 
by definition giving 2 J. This problem is set up in Ref. [l3l 
and consists of constructing wave functions ifjoi 2 ) an d 
4>o(—z) localized in the left and right wells, respectively. 
It is then shown in a straightforward way that the energy 
splitting in this double well is given by 



J 



-^(O)Vo(O). 



(10) 



To evaluate this expression, it is crucial to distinguish 
cases of a smooth and sharp barriers. The case of a 
smooth barrier is solved in Ref. 13 using the WKB ap- 
proximation and leads to the expression for J propor- 
tional to the natural frequency of the oscillations (at- 



tempt frequency) in each of the wells as well as the dimen- 
sionlcss coefficient of penetration through the barrier. In 
contrast, motivated by the connection to the J J problem 
we are instead interested in a sharp barrier (relative to 
the size L of the wells) of width d, located between points 
z = —d/2 and z = d/2 (z is the axis along the leads and 
perpendicular to the barrier), where the potential energy 
exceeds the kinetic energy by the amount Uq. For this 
setup the WKB is clearly inapplicable. Under such con- 
ditions, the wave function in the wells is approximately 



(11) 



where k is the wave vector in the wells, while the wave 
function under the barrier is 

U*) « -U^^' -± < z < 0, (12) 
v L ^ 



with length I = fij \ / 2m a Uo. This then gives 14 

h 2 



J 



-d/e 



JL 



(13) 



Another more direct way of producing Eq. (|13p is to 
note that J is proportional to the matrix element of the 
Hamiltonian (e.g., the kinetic energy) between the left 
and right wavefunctions, normalized in each of the wells 
and penetrating a distance £ into the barrier. Clearly 
then J w J^ L dz^-ip^ip'^ giving result (jT3j>, and as ad- 
vertized vanishes for macroscopically large wells. 

Consequently, for any realistic energy difference 2h be- 
tween the two wells, J in ([8]) can be neglected, reducing 
the Rabi oscillation frequency for such large wells to 



lor = 2h/H, 



(14) 



and allowing the identification of the imbalance h in the 
double-well BEC with the chemical potential difference 
(voltage) eV, ©. 

We now turn to the analysis of the effects of interac- 
tions on our conclusion above by considering the Hamil- 
tonian ([?]) with gi 7^ 0. This is most conveniently done in 
the (polar) density-phase representation, by introducing 



and 



N 

N L = -(l + m), 



(15) 



(16) 



N 



(l-m), cj>! 



e 



where — 1 < m < 1. The total number of particles N is 
conserved (hence H is independent of 9), and the part of 
the Hamiltonian describing the evolution of m becomes 



H 



H =jf = -Jv 7 ! - m 2 cos(20) + hm + ^m 2 , (17) 
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where A = (go — <7i)/2. The equations of motion are thus 



—h — Xm — J - 
dU 



to , du , . 

jC os(2<£) = ^, (18) 



vl — m 



hrn = — — = 2J\J 1 — to 2 sin(20), 



(19) 



with (fl8|) defining the effective chemical potential differ- 
ence Sfi = ^l—\ 1 r- Although it is not fixed in this canon- 
ical ensemble, in the thermodynamic limit of interest its 
constant (in time) part can be meaningfully associated 
with the grand-canonical chemical potential different of 
the JJ system. Above equations provide a complete solu- 
tion to the problem of two couple interacting condensates 
and have been analyzed in Ref. [f| Although it is possible 
to solve these equations analytically, the general solution 
is not very informative. 

Here we focus on the physical limit of large reservoirs 
discussed above, in which the coupling J is vanishingly 
small. In this limit we solve Eqs. (|18[) perturbatively in 
powers of J/h <C 1. At zeroth order in J/h, 



h(f)Q = —h — XrriQ, 
to = 0, 



leading to the solution 



to 



h + Ato,q 



-ht. 



const. 



Identifying atom current I with (Nr — Nl)/2 
we obtain, with the help of Eq. (fi"9| , 



JNy/l 



2(h + Xm )t 



(20) 
(21) 



(22) 
(23) 

-Nm/2, 
(24) 



that is clearly equivalent to with the identification 
© and I = JNy/l - mj/h. 

From this last identification of Iq we can further 
note that the result Eq. (fT3|) is compatible with the 
Ambegaokar-Baratoff expression for the critical current 
in a Josephson junction, namely, that 



N 
L£ 



pJl. 



(25) 



with p = N/L 3 the atom density. From this we observe 
that the critical current of the double- well BEC scales as 
the area L 2 of the "junction" (double- well barrier), just 
like the Ambegaokar-Baratoff expression for the super- 
conducting JJ, Iq B = 7rAG„/(2e), where the junction 
area enters through the number of conduction channels 
in G n . This further supports our finding that J, as de- 
fined in ((J) vanishes with the inverse length of the well 
"leads". 



Let us see whether any realistic Josephson junctions in- 
deed obey the condition J <§; /i. For a typical Josephson 
junction, d ~ lOOnm [l5[. I must itself be of the order 
of d, otherwise the coupling J will be even further expo- 
nentially suppressed. Taking L > lOOnm and h ~ leV, 
we find 



J 
h 



m a £Lh 



< 10" 



(26) 



where an electron mass was used for m a . Roughly, the 
above factor of 10 -5 arises due to two factors of 10 3 of 
I and L (large characteristic length scales) relative to 
the Bohr radius that corresponds to an eV energy scale. 
Thus indeed in a realistic Josephson junction J is always 
much smaller than h, well justifying above approxima- 
tion, even if the system were to be closed. Of course, 
as mentioned in the introduction, a conventional super- 
conducting Josephson junction is an open system, that 
is a part of a macroscopic circuit and is therefore effec- 
tively characterized by an infinite L. Furthermore, with 
the circuit driven by a fixed voltage source the Josephson 
frequency expression ujj, |(5J) is effectively exact; any in- 
accuracy in Josephson voltage-frequency relation quoted 
above is associated with the uncertainty of the current 
knowledge of Planck's constant (one part in 10 7 )[l6l|. 

In contrast a typical double-well BEC trap potential 
is expected to be roughly characterized by a single (the 
same order of magnitude) length and energy scale, with 
i ~ L and Uo ~ h, leading to J and h, that are compa- 
rable and both a tiny fraction of an electron volt. Con- 
sequently, we expect such a system to display a signifi- 
cant and tunable deviation from the Josephson frequency, 
Eq. (j2"2")l , obtained by neglecting J/h corrections. 

Even though as discussed above in a double-well BEC 
system we generically expect J ~ h, we observe that the 
experimentally studied double- well BEC0 is character- 
ized by J/h w 15 sec" 1 , and 1400 < h/h < 5700 sec" 1 
(here and throughout h denotes the chemical potential 
imbalance energy, not the Planck's constant 2irh), and 
thus corresponds to 0.003 < J/h < O.Oll.pj This small 
value of J/h characterizing these experiments explains 
why the measured oscillation frequency (u> ac .j) as a func- 
tion of the chemical potential difference (h) is observed 
to be linear in Ref. d. In order to detect a deviation from 
this linear behavior these measurements need to be ex- 
tended down to h/h < 15 sec -1 , or done on a system in 
which the value of J is increased by e.g., making wells 
smaller. 

With this in mind, it is useful to compute the lowest 
order correction. This can be straightforwardly done by 
evaluating the solutions <j)(t), m(t) to Eq. (fT8")) systemati- 
cally to nth order in J/h by iterating the equations, with 
the nth-order solution <j> n (t) , m„ (t) on the left-hand side 
and approximating the right-hand side by the n — 1st- 
order solution, <fi n _i(t),m n -i(t). 



4 



To first-order the equations become: 



b± = -h-J 



to 



:COS(20 O ), 



toi = 2 Jy 1 — TOq sin(2</>o 
leading to the solution 

J toq 



-ht 



2h y/l - ml 



sin(2M), 



J 



TOi(i) = mo + ^- \/l — to 2 cos(2fti). 
h 



The solution to second-order, (J/h) 2 is obtained by 
using above expressions, </>i(i), mi(t) in the right-hand 
(27) side of the exact equations (fl~8|) : 



(28) 

(29) 
(30) 



-h-J 



mi 



■. cos(2(/>i), 



TO 2 



2JW1 - TOf sin(2<^i) 



(31) 
(32) 



Focusing on <p(t) and integrating the first equation we 
find the time-independent part of <f> 



-h - J- 



toq + 6m(t) 



: COS 



2ht + 'i- 



m 



Jtoq 



: cos(2ht) 



J 



h 1 - to; 



■sin 2 (2/rf) 



: sin(2M) 



/i 1 - to; 



■cos 2 (2/if), 



J 2 



J 



2/i 1 - mi 2hl-mi 



4«-\//' 2 
2/i 



J 2 , 



(33) 
(34) 
(35) 
(36) 



which by definition ([18]) is the effective chemical potential 
difference. In above we have neglected the higher order 
hamonics that are also always generated at nonzero J/h. 
This then gives 



TO 



t^h 2 + J 2 , 



2J\ 1 — to 2 , sin 



2t\ h 2 + J 2 



(37) 
(38) 



We therefore obtain the (fundamental) ac Josephson fre- 
quency of current oscillations (defined by to) to be given 
by 



I j2 
h+^r 

2h 



h 2 + J 2 , 



(39) 



with h = h + Atoo. 

Thus, as advertized, for macroscopically occupied wells 
(the only limit in which a reasonable connection to a 



Josephson junction can be made), such that J/h <C 1, 
an imbalanced double-well BEC system is indeed a good 
model for the ac Josephson effect, exhibiting current os- 
cillations with frequency that is nearly independent of J 
and grows linearly with imbalance h — h+Xmo- However, 
for a smaller double-well BEC (more typical experimen- 
tally) the fundamental oscillation frequency is expected 
to exhibit (J/h) 2 deviations from the linear dependence 
on h of the form given in (|39|) . 
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